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The quantum phase transition between paramagnetic and antiferromagnetic phases of the Kondo 
lattice model with Ising anisotropy in the intersite exchange is studied within extended dynamical 
mean-field theory. Nonperturbative numerical solutions at zero temperature point to a continuous 
transition for both two- and three-dimensional magnetism. In the former case, the transition is 
associated with critical local physics, characterized by a vanishing Kondo scale and by an anomalous 
exponent in the dynamics close in value to that measured in heavy- fermion CeCu5.9Auo.1- 



PACS numbers: 71.10.Hf, 71.27. +a, 75.20.Hr, 05.10.Cc 

Heavy-fermion materials close to a zero-temperature 
antiferromagnetic (AFM) instability manifest a rich vari- 
ety of non- Fermi- liquid behaviors 1] . By tuning a control 
parameter (doping, pressure, or magnetic field), the Neel 
temperature can be suppressed to T = 0. The result- 
ing quantum phase transition (QPT) separates an AFM 
metal from a paramagnetic (PM) metal in which local 
moments are Kondo-screened by conduction electrons. In 
the standard spin-density wave (SDW) theory [2| for such 
QPTs, the important critical modes are long-wavelength 
order-parameter fluctuations; non-Fermi-liquid physics 
arises from scattering of quasiparticles by SDWs. How- 
ever, data for a number of heavy fcrmions are at odds 
with this picture. In particular, neutron-scattering for 
CeCus.gAuo.i [3j,|4( points to a strongly interacting quan- 
tum critical point involving novel local physics [5(. 

Many of the deviations from the SDW theory can be 
explained if criticality in the magnetic ordering renders 
Kondo physics simultaneously critical, producing a "lo- 
cally critical QPT" @, 0, H- 

Here, work has focused 
on the Kondo lattice model, mapped within extended 
dynamical mean- field theory (EDMFT) @ onto a self- 
consistently determined Bose-Fermi Kondo (BFK) model 
for a magnetic impurity coupled both to a fermionic band 
and to a dissipative bosonic bath representing nonlo- 
cal fluctuations arising from Ruderman-Kittel-Kasuya- 
Yoshida (RKKY) interactions between local moments. 

The central question we set out to answer is whether 
the EDMFT captures a continuous QPT. Such a transi- 
tion is expected to be of SDW type for three-dimensional 
(3D) magnetic spin fluctuations, with local moments that 
are Kondo-screened below a finite low-temperature scale 
Tk and hence play no role at the QPT. By contrast, if the 
fluctuations are two-dimensional (2D) — for which there 
is some evidence in CeCus.gAuo.i [1Q|] — the scale is 
claimed 0] to vanish precisely at the QPT with coexist- 
ing long-wavelength and local critical modes. To confirm 
the scenarios described above, it must be demonstrated 
that the PM and AFM phases terminate at a common 
value of the control parameter. Numerical studies em- 
ploying quantum Monte Carlo impurity solvers have ad- 
dressed the particular case of Ising anisotropy believed 
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The character 
11] , in part due 



to be relevant to CeCu,' 
of the QPT remains contentious [7|, 
to the limitations of quantum Monte Carlo methods in 
accessing the lowest temperature scales. Solutions of the 
Kondo lattice model, found via a mapping of the effective 
BFK model onto the spin-boson model, have pointed to 
a continuous, locally critical QPT for the 2D case jjjj]. By 
contrast, studies of the related periodic Anderson model 
have predicted a strongly first-order QPT, casting doubt 
on the ability of EDMFT to capture the key physics 11 1 . 

This Letter presents zero-temperature EDMFT solu- 
tions, obtained using a numerical renormalization-group 
(NRG) technique [12| to solve the underlying BFK 
model. Results for both 2D and 3D spin fluctuations 
point to a continuous QPT, at which the magnetic order 
parameter vanishes continuously as the static suscepti- 
bility at the AFM ordering wave vector diverges, and 
the terminus of PM solutions coincides with the onset of 
magnetic ordering. For the topical 2D case, a concomi- 
tant divergence of the static local susceptibility signals 
that Kondo physics is driven critical at the QPT, and 
we are able to extract an anomalous exponent in good 
agreement with experiments on CeCus.gAuo.i [J]. These 
results provide important support for the notion of lo- 
cal quantum criticality @> 0> @] ' both as a description 
of specific heavy-fermion systems and, more generally, as 
a paradigm for novel phase transitions in other areas of 
physics. 

Model. — We focus on the Kondo lattice model with 
Ising RKKY anisotropy, specified by the Hamiltonian 
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A spin-i local moment Si at site i is coupled to the on- 
site conduction-electron spin s c i via an exchange J > 0, 
favoring quenching of local moments. The competing 
tendency towards magnetism enters through the RKKY 
interaction 7y between the z components of the local- 
ized spins on different sites i and j. The tight-binding 
parameters tij determine the conduction-band dispersion 
efc and hence the density of states p(e) = J2k — £fc )- 
EDMFT formulation— We study the Kondo lattice flU 
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using EDMFT 0] , which amounts to considering a single 
lattice site, described by the impurity Hamiltonian 

#BFK = ^ C ka C ka + wM>q + JS-S C 



Sz J2 9i U>q + <f>-q) + hioA, (2) 



k.a 



at which a single spin S interacts with the on-site spin s c 
of a conduction band and, via S z , with both a dissipative 
bosonic bath (representing magnetic fluctuations due to 
spins at all other lattice sites) and a local magnetic field 
h\ oc . The coupling of the impurity to the nonlocal degrees 
of freedom is fully specified by Jp(e) and the bosonic 
spectral function B(uj) = -k g q 5(u> — u) q ). 

Within EDMFT, the dynamical spin susceptibility is 
written x(q,a>) = [M(u>) + Iq}^ 1 in terms of a "spin 
self-energy" M(u>) and I q , the spatial Fourier transform 
of lij. The key approximation in EDMFT is neglect 
of the momentum dependence of both M and the self- 
energy E CT entering the conduction-electron Green's func- 
tion G a (k,e) = [e — £fc — E^e)] -1 . M and are 
determined self-consistently by requiring that the wave- 
vector average of the lattice correlation function equals 
the corresponding local correlation function of the im- 
purity problem ((2]). Thus, the local susceptibility must 
satisfy 



Xioc(^) = i dte 1 
Jo 



t ([S z (t),S z (0)}) = 



Me) de 
M (w) + e 



(3) 



ql being the RKKY density of states. 



The Weiss field Xo (w) — M(uS) — Xioc( w ) determines 
the bosonic bath density of states through Imx ( ^ 1 (ti>) = 
sgnHi?(M). Analogously, G^(e) = G^Je) - E a (e) 
determines the band density of states p{e). However, to 
study critical properties, it is not necessary to enforce 
self-consistency on G a (k, e). Following [7], we instead 
take p(e) = po 0(D— |e|) and work with fixed poJ, taking 
D = 1 as the energy unit. This choice of a feature- 
less conduction band avoids double counting of RKKY 
interactions, discussed 0] as the source of discrepancies 
between previous T > studies 0, @, 11 1 ■ 

For PM solutions, h\ oc = 0, and we follow solutions 
upon increasing the control parameter S = I /T^-', here, 
— / = Iq is the most negative value of I q , found at the 
AFM ordering wave vector q = Q, and = Xioc( w = 
0; 6 = 0) is the bare Kondo scale of the pure-fermionic 
impurity problem. Instability of the PM phase is signaled 
by a divergence of x(Q> 0) f° r some S — <5 Ci p. For AFM 
solutions, h\ oc is related self-consistently to a nonzero 
staggered magnetization /tt-afm = (S z ) by 



hlac = - [I - Xo 1 = 0)] TOAFM • 



(4) 
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We follow AFM solutions from large 5 down to S — <5 c .a, 
below which the only solution is ttiafm = 0. 



1: (Color online) Converged values of y = [/xioc(0)] 1 
vsx = [/x(<5,0)] _1 . In 3D (squares), Eq. © forces y = |(1 + 
x) + | \/ x(x + 2) (dashed line) to approach \ as x vanishes. 
In 2D (circles), both quantities vanish together according to 
y — 2/ ln(l+2/a;) (solid line), signaling critical Kondo physics. 



Examination of the self-consistency condition Q at 
u> = reveals that the nature of the QPT depends cru- 
cially on the dimensionality of the spin fluctuations: In 
3D, pi(e) has a square-root onset at its lower edge e = —I. 
We take Pl (e) = 2VP - e 2 /{nl 2 ) 0(1 - |e|), for which 
Eq. (O gives xioc(^) = 4%o 1 ( LJ )/I 2 i an d the local static 
susceptibility xi oc (0) remains finite when the peak sus- 
ceptibility x{Q,0) diverges (dashed line in Fig. Q]). For 
2D magnetism, p\(e) instead has a jump onset, carica- 
tured as pi(e) = (2I)- 1 e{I~ \e\), for which Eq. © yields 
Xioc(cj) = (2/) _1 In [1 + 2Ix{Q, w)]- In this case, any di- 
vergence of x(Q, 0) is necessarily accompanied by a di- 
vergence (albeit weaker) of Xioc(0) (solid line in Fig. |T|) . 

To determine where on the parametric curves shown 
in Fig. Q] the EDMFT solution lies for a given S, it 
is necessary to solve the full set of EDMFT equations 
at all frequencies. The key issue is whether, for 2D 
magnetism, EDMFT yields a locally critical QPT, i.e., 
whether S c ^a — <5c,p = <5 C with Xioc^ — 0; 5 — S c ) = 0. 

Solution method. — We solve Eq. (|2|) using an extension 
(l2| of the NRG [l]| . After logarithmic discretization of 
the energy axis (|e|, to = DA~ n for n = 0, 1, 2, . . . with 
A > 1), Hbfk is recast in terms of a fermionic and a 
bosonic chain, each coupled to the impurity at the first 
site only. An exponential decay of tight-binding coeffi- 
cients along each chain allows iterative diagonalization 
of chains of increasing length. The Fock space must be 
truncated, allowing a maximum of Nb bosons per site, 
and retaining only the lowest N s many-body eigenstates 
from one iteration to construct the basis for the next. 

The NRG method gives an excellent account of the 
universal critical properties of the BFK model, which 
converge rapidly with increasing Nb and N s , and are 
insensitive to A [T3|. However, obtaining a satisfac- 
tory description of the Kondo lattice places more strin- 
gent demands on the impurity solver. The NRG gives 
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the imaginary part of Xioc(w) as a set of delta func- 
tions that are broadened to recover a continuous Xioc( w ) 
[lH; the real part x( oc (w) follows by Hilbert transforma- 
tion. When h loc — ► 0, Xloc(0) so obtained should co- 
incide with —(S z }/h\ oc , but in practice there is a mis- 
match between these two values. If the self-consistently 
determined itiafm, and hence h\ oc , vanish continuously 
at 8 c ,a, then by rewriting Eq. |4]) in the equivalent form 
^ioc/"iAFM + Xioc(0) = X _1 (QjO), one sees that any such 
mismatch prevents divergence of x(Q>0)j which in turn 
leads to coexistence of PM and AFM solutions. It is this 
previously neglected mismatch that underlies the uncon- 
ventional scenarios presented in 15, la ] , and not trunca- 
tion error or a failure of EDMFT as there suggested. 

We ensure internal consistency of the NRG by apply- 
ing a multiplicative correction 1 + c to xioc(<^, ^loc), such 
that xioc(0, rj) — ~(S z )/i], where ij « 10~ 6 is the small- 
est field used in our calculations. For A = 3, c » 0.15 for 
the BFK model, very similar to the value for the pure- 
fermionic Kondo model. Within the EDMFT, the static 
lattice susceptibility x(Q, 0) now diverges by construc- 
tion if h\ oc — > r\ self-consistently, i.e., if the order param- 
eter vanishes continuously — a scenario that we emphasize 
is in no way imposed by the correction scheme. 

Self-consistent EDMFT solutions are obtained as fol- 
lows: For given /, Eq. ((2|) is solved using the NRG with a 
trial B(lo) and an initial field h\ oc > rj [h\ oc = rj\ for AFM 
[PM] solutions. At this point the numerical mismatch is 
corrected. For AFM solutions, finding c requires an ad- 
ditional NRG calculation performed at h\ oc — r\. Then, 
B{tu) [and for AFM solutions, h\ oc > rj\ is updated for 
use in the next loop via Eq. Q [and Eq. (J4])]. Typi- 
cally 10-50 EDMFT loops suffice to converge B(oj) to 
within 0.0001% for all ui, though convergence is consid- 
erably slower around the critical couplings. The solution 
reached is independent of the details of the trial B(ui) 
and, for AFM solutions, of the initial h\ oc . 

Results are presented for A = 3, Nb — 8, and N s — 300 
(requiring up to 1 hour per EDMFT loop on a 2.2-GHz 
AMD Opteron CPU). All data shown are for p a J = 0.2, 
corresponding to a bare Kondo scale w 0.014, but 
we find near-perfect scaling with Tjj for T^- < 0.1. Wc 
have established in representative cases that the numerics 
are converged with respect to Nb and N s , and that the 
essential physical picture is independent of A. 

Results for 3D magnetism. — In this case, EDMFT has 
been supposed (but not explicitly shown) to exhibit a 
continuous QPT of the SDW type p. This picture 
is confirmed in Fig. Eta), which shows the variation of 
three T = static quantities with S = I/T^. The 
inverse peak susceptibility X^HQjO) of the PM solu- 
tions decreases with increasing S, and vanishes linearly 
at 6 CjP = 1.33(2) [see Fig. E^a)], beyond which no PM 
solutions are found. The local susceptibility xioc(0) ap- 
proaches Tj^xi oc (0) = 2/<5 c .p at <5 c .p, as follows from the 
self-consistency (see Fig. [T]). For large 5, AFM solutions 
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FIG. 2: (Color online) Peak and local static susceptibilities 
(in units of T^) of PM and AFM solutions, along with the 
order parameter toafm, for (a) 3D and (b) 2D magnetic 
fluctuations. AFM solutions for <5 = I /Tf { < <5 c ,a (where 
?tiafm = self-consistently) coincide with PM solutions (with 
7TIafm = 0). Lines are guides to the eye. 



converge with nonzero h\ oc and toafm- Our data are 
consistent with a continuous vanishing of toafm upon 
decreasing S and, provided the internal consistency of 
the NRG is enforced (see Solution method), a concomi- 
tant divergence of x(Qj0)- Extrapolation of the low- 
est x _1 (Qj0) values to zero yields (5 Cj a = 1.31(2). For 
S < S c _a, AFM and PM solutions coincide. 

The large (around 30% of <5 c ,a) region of coexisting 
PM and AFM solutions reported in [l5( is sharply sup- 
pressed by enforcing the internal consistency of the NRG. 
The reduced coexistence range of 0-5% likely stems from 
residual errors (due to discretization, truncation, Hilbert 
transformation, numerical rounding, etc.) that will tend 
to destabilize the fine balance inherent to any continuous 
transition. We conclude that, within numerical accuracy, 
the EDMFT description of the QPT for 3D magnetism 
is continuous and of conventional SDW type Q . 

Results for 2D magnetism. — Static quantities for the 
2D case are shown in Fig. [2fb). Convergence slows 
markedly in the vicinity of the transition due to the 
logarithmic form of the self-consistency, making it ex- 
ceedingly difficult to obtain solutions having very small 
X _1 (Q,0). PM-phase calculations also require a careful 
choice of starting parameters to avoid flowing towards an 



unphysical solution 14]. Nonetheless, linear extrapola- 
tion of X 1 (Q: 0) to zero for each phase [Fig. EJa)] yields 
5 C , P = 1.40(2) close to S c ,a = 1.34(2). The ground-state 
energy (not shown) provides another handle on the po- 
sition of the transition. Extrapolations of the PM and 
AFM energies indicate a crossing at S c — 1.36(1), consis- 
tent with <5 C A deduced from x _1 (Q ; 0)- 

The slower convergence in 2D has prevented us from 
approaching the transition as closely as in 3D, or from 
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FIG. 3: (Color online) (a) Extrapolation to zero of PM 
X _1 (Q,0) (squares) and AFM x _1 (Q,0)/20 (circles) for 3D 
and 2D magnetism, (b) Real part of Xioc(^) for five 8 values 
approaching 8 C in the 2D AFM phase. The dashed line is a fit 
to the logarithmically singular form that develops for 8 — > S c . 



temperature solutions of the Kondo lattice model, of 
great present interest in connection with heavy-fermion 
quantum criticality. For magnetic fluctuations both in 
two and three dimensions, our results point to a contin- 
uous transition. In the 2D case, critical local-moment 
fluctuations are observed with an anomalous exponent 
a ~ 0.8 in the dynamics that is in good agreement with 
the experimentally determined value for CeCus.gAuo.i. 
This provides significant new evidence for local quantum 
criticality in strongly correlated systems. 

We thank R. Bulla, P. Coleman, S. Kirchncr, P. Ku- 
mar, Q. Si, P. Sun, P. Wolfle, and J.-X. Zhu for stimulat- 
ing discussions. We acknowledge use of the University of 
Florida High-Performance Computing Center and thank 
C. Taylor for technical support. This work was supported 
in part by NSF Grant DMR-0312939. 



converging AFM solutions below ttiafm ~ 0.14 (after 
more than 300 EDMFT loops). However, the narrow co- 
existence region (1-8% of S c ,a) is comparable with that 
in 3D, and to within our numerical accuracy is consis- 
tent with a continuous transition. Such a QPT is nec- 
essarily of the locally critical type @, with xioc(0) and 
x{Qi 0) diverging together at S c (see Fig. [1]), where the 
self-consistently determined BFK model lies at its critical 
point [17]. We cannot rule out a very weakly first-order 
transition either in 2D or 3D, but physically this scenario 
is virtually indistinguishable from the continuous case. 

Support for critical local physics is provided by Xiod 10 ) 
[Fig. GH[b)], which develops a logarithmically singular 
form H T°.x(ocH ~ (a/2S c ) In \oj\~ 1 in the vicinity of 
the critical coupling for T£ -C \u>\ -C 1%, where the ef- 
fective Kondo scale Tj£- vanishes logarithmically slowly as 
S -► S c . Wc find a/28 c = 0.29(1) or a = 0.78(4), consis- 
tent with the value measured in CeCus.gAuo.i (a — 0.75) 
and with that obtained in one T > study No such 
anomalous exponent was found in connection with the 
strong first-order behavior reported in [H| . 

A separate EDMFT study of the Kondo lattice has 
reached similar conclusions about the nature of the QPT. 
In [l^ . the fermionic degrees of freedom are eliminated 
by mapping the BFK model onto a spin-boson model. 
While this mapping introduces another level of approx- 
imation and rules out direct comparison of the ground- 
state energy in the PM and AFM phases, the spin-boson 
model can be solved via a purely bosonic NRG method 
requiring less computational effort per EDMFT iteration, 
allowing investigation closer to the critical points where 
convergence slows markedly. That two complementary 
approaches yield essentially the same physical picture is 
an important confirmation of the key results. 

Summary. — We have obtained nonperturbative, zero- 



[1] 
[2] 

[3. 



[6] 

[7] 

[8 
[9 

[10 
[11 

[12 

[13 

[14; 

[15 
[16 

[17; 

[18 



G. R. Stewart, Rev. Mod. Phys. 73, 797 (2001). 

J. Hertz, Phys. Rev. B 14, 1165 (1976); A. J. Millis, 
Phys. Rev. B 48, 7183 (1993). 

H. v. Lohneysen et ai, Phys. Rev. Lett. 72, 3262 (1994); 
A. Schroder et ai, ibid. 80, 5623 (1998); O. Stockert et 
ai, ibid. 80, 5627 (1998). 

A. Schroeder et ai, Nature (London) 407, 351 (2000). 
P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 
Phys. Condens. Matter 13, R723 (2001); P. Coleman and 

C. Pepin, Physica B 312, 383 (2002); Acta Phys. Pol. B 
34, 691 (2003). 

Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Na- 
ture (London) 413, 804 (2001); Phys. Rev. B 68, 115103 
(2003). 

D. R. Grempel and Q. Si, Phys. Rev. Lett. 91, 026401 
(2003); J.-X. Zhu, D. R. Grempel, and Q. Si, Phys. Rev. 
Lett. 91, 156404 (2003). 

Q. Si, J.-X. Zhu, and D. R. Grempel, J. Phys. Condens. 
Matter 17, R1025 (2005). 

Q. Si and J. L. Smith, Phys. Rev. Lett. 77, 3391 (1996); 

J. L. Smith and Q. Si, Phys. Rev. B 61, 5184 (2000); R. 

Chitra and G. Kotliar, Phys. Rev. Lett. 84, 3678 (2000). 

A. Rosch et ai, Phys. Rev. Lett. 79, 159 (1997). 

P. Sun and G. Kotliar, Phys. Rev. Lett. 91, 037209 

(2003); Phys. Rev. B 71, 245104 (2005). 

M. T. Glossop and K. Ingersent, Phys. Rev. Lett. 95, 

067202 (2005); Phys. Rev. B 75, 104410 (2007). 

R. Bulla, T. Co sti, and T. Pruschke, 

|arXiv:cond-mat/0701105| 

K. Haule, A. Rosch, J. Kroha, and P. Wolfle, Phys. Rev. 
Lett. 89, 236402 (2002). 



M. 



Glossop 



and 



K. 



arXiv:cond-mat/0607566vl 
Kirchner 



R. Bulla, 



Ingersent, 
and Q. Si, 



j.-x. zhu, s: 

arXiv:cond-mat/0607567vl. 

Refs. Q describe, for 8 > 8 C and T > 0.01T&, "nominally 
PM" solutions belonging to the localized phase of the 
effective BFK model. We find no such solutions at T — 0. 
J.-X. Zhu, S. Kirchner, R. Bulla, and Q. Si, following 
Letter, Phys. Rev. Lett. 99 227204 (2007). 



